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The location of the critical end point of QCD has been determined in previous studies of Nf = 2 + 1 
and A/ = 2 + l + l dynamical quark flavors using a (truncated) set of Dyson-Schwinger equations 
for the quark and gluon propagators of Landau-gauge QCD. A source for systematic errors in these 
calculations has been the omission of terms in the quark-gluon interaction that can be parametrized 
in terms of baryonic degrees of freedom. These have a potentially large dependence on chemical 
potential and therefore may affect the location of the critical end point. In this exploratory study 
we estimate the effects of these contributions, both in the vacuum and at finite temperature and 
chemical potential. We find only a small influence of baryonic contributions on the location of the 
critical end point. We estimate the robustness of this result by parameterizing further dependencies 
on chemical potential. 


I. INTRODUCTION 

Heavy ion collision experiments at RHIC/BNL and the 
future CBM/FAIR facility probe the phase structure of 
QCD at finite chemical potential. One of the major goals 
of these experiments is the study of the existence, the 
location and the properties of a critical end point (CEP), 
where the chiral crossover at small chemical potential 
turns into a first-order transition. 

From a theoretical point of view it remains unclear at 
the moment whether this is indeed the case. Lattice cal¬ 
culations firmly established the crossover behavior at zero 
chemical potential, see e.g. [1, 2] and references therein. 
At finite chemical potential, lattice calculations are ham¬ 
pered by the notorious fermion sign problem. Although 
various extrapolation methods agree with each other at 
rather small chemical potential [3-11], for regions in the 
(T,ff) plane with pb/T > 2 uncertainties accumulate 
rapidly. Therefore, despite many efforts the basic prop¬ 
erties of the phase diagram of QCD are not yet settled, 
see e.g. [12, 13] and references therein. Thus other theo¬ 
retical methods are mandatory to complement the lattice 
calculations. 

In a sequence of previous papers [14-16] the loca¬ 
tion of the CEP tor Nf = 2, Nf = 2 + 1 and Nf = 
2 + 1 + 1 quark flavors has been determined using the 
(truncated) Dyson-Schwinger equations (DSEs) for the 
quark and gluon propagators of QCD. A particular fo¬ 
cus in these studies has been the inclusion of back- 
reaction effects of the quarks onto the Yang-Mills sec¬ 
tor, which allowed to go beyond simple modeling of 
the gluon part either within the DSE approach [17, 18] 
or in chiral models like the Nambu-Jona-Lasinio (NJL) 
model [19], its Polyakov-loop extended versions [20-22] 
and the Polyakov-loop extended quark-meson (PQM) 
model [23-25]. At zero chemical potential the inclusion 
of the quark back-reaction on the gluons produced the 
correct temperature behavior of the quark condensate 
and led to predictions for the magnitude of unquench¬ 
ing effects in the gluon propagator [14] , which have been 
verified by subsequent lattice calculations [26]. 


In this approach a CEP has been found at rather large 
quark chemical potential (T°,/rp = (115,168) MeV. 
Since this result relies on a truncation of the quark-gluon 
interaction which is still far from complete, it is an impor¬ 
tant task to quantify its systematic error. In particular, 
effects due to a nonzero chemical potential that cannot be 
tested by comparison with lattice calculations at = 0 
may provide for sizable quantitative corrections. In this 
work we focus on a particular class of such corrections, 
namely vertex corrections that can be parametrized in 
terms of (off-shell) baryons. In general, baryonic back- 
reaction effects onto the quark propagator provide a di¬ 
rect mechanism how the quark condensate may be influ¬ 
enced by changes in the baryon’s wave functions such as 
the one inflicted e.g. by the nuclear liquid-gas transition 
at very small temperatures [12]. These back-reaction ef¬ 
fects, however, may very well decrease in size for growing 
temperatures and it is an open question whether they are 
still important in the region of the QCD phase diagram 
where the putative CEP for the chiral phase transition 
is located. In a two-color version of QCD this influence 
has been studied in Refs. [27, 28] and found to be crucial 
to an extent that not only the location but even the very 
existence of a CEP is affected. Whether this is also the 
case in the SU{3) theory is an open question that needs 
to be addressed. We regard the study reported in this 
work as a first step in this direction. 

The paper is organized as follows. In Sec. II we re¬ 
view our truncation scheme of the DSEs for the quark 
and gluon propagators. We explain the details of the 
corresponding equation for the quark-gluon vertex and 
identify the diagrams that can be parametrized in terms 
of (non-elementary, i.e. composite) hadronic degrees of 
freedom. We specify an approximation scheme for these 
terms that is suitable for an exploratory calculation of its 
effects on dynamical chiral symmetry breaking. For sim¬ 
plicity we restrict ourselves to the case of QCD {Nc = 3) 
with two degenerate fermion flavors, Nf = 2, and note 
that a generalization to the Nf = 2 + 1 case is straight¬ 
forward but very expensive in terms of CPU time. In 
Sec. Ill we present our results. We first discuss the effects 
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of the resulting baryon loop on the quark propagator in 
the vacuum and at hnite temperature but zero chemical 
potential. We then present an estimate for the size of the 
effects that may be expected for the CEP. We conclude 
in Sec. IV. 


II. DYSON-SCHWINGER EQUATIONS 
A. DSEs for the propagators 


In order to accommodate the notation already for its 
intended purpose later in this work, we specify the quark 
and gluon propagators at finite temperature T and quark 
chemical potential Hq and indicate the limits T —>■ 0 and 
/j,g —0 where appropriate. The bare quark propagator 
is then given by Sq^{p) = ip ■ j Z2 + ioJnlA ^2 + ^2 w, 
with wave function renormalisation Z2 and bare quark 
mass TO. The dressed inverse quark propagator S~^ and 
the Landau-gauge gluon propagator are given by 


S 1 (p) = fp • 7 A{p) A iCbnli C{p) + B{p), 

d,Ap) = + pUp)^ 


( 1 ) 


with momentum p = {ujn,p)- The Matsubara frequencies 
are a;„ = 7rT(2n -|- 1) for fermions and a;„ = 7rT2n for 
bosons, and we use the abbreviation = a;„ -I- ipq. All 
dressing functions implicitly depend on temperature and 
chemical potential. The projectors are transverse 
(T) and longitudinal (L) with respect to the heat bath 
and given by 

P^v = (1 - (1 - ^ 

pL _ p _ pT 

where = 5^^, — p^Pv/p^ is the covariant transverse 
projector. The limit of zero chemical potential is straight¬ 
forward; in the additional zero-temperature limit the 
momentum p reduces to its usual 0(4)—symmetric Eu¬ 
clidean form and the wave functions of the quark prop¬ 
agator become degenerate, i.e. A(p) = C'(p). Further¬ 
more, the transverse (magnetic) and longitudinal (elec¬ 
tric) dressing functions of the gluon approach the same 
limit in the vacuum, i.e., Zt{p) = ■Z^l(p) = Z{p). In 
the medium there exists also a fourth contribution to the 
inverse quark propagator, which vanishes in the vacuum 
and is proportional to uinjAP ■ 7- Due to its negligible 
contribution also at higher temperatures and chemical 
potential we do not consider it throughout this work. 

The DSE for the quark propagator is shown diagram- 
matically in Fig. 1. The pieces that need to be de¬ 
termined in order to allow for a self-consistent solution 
of this equation are the fully dressed gluon propagator 
and quark-gluon vertex. Model calculations [17, 18] of¬ 
ten use simple ansdtze for the gluon propagator that do 



FIG. 1. The DSE for the quark propagator (top panel) and 
the truncated gluon DSE for Nf = 2 QCD (bottom panel). 
Large blobs denote dressed propagators and vertices, and the 
hatched circle represents the quenched (lattice) propagator. 


not take into account the proper temperature and fla¬ 
vor dependence of the gluon self-energy. We prefer to 
include these important effects by taking the Yang-Mills 
sector of QCD into account and calculating the back- 
reaction of the quarks onto the gluon explicitly. This 
framework has been gradually evolved from the quenched 
case, Nf = 0 [32, 33], to two-flavor QCD [14, 34] and re¬ 
cently to Nf = 2 + 1 and TV/ = 2 -|- 1 -|- 1 [14, 16]. 

Such an approach has two distinct advantages over 
simple modeling. On the one hand it allows us to trace 
the effects of quark masses and flavors as exposed in the 
Columbia plot. On the other hand, it serves to take into 
account the effects of chemical potential on the gluon 
propagator explicitly, thereby rendering results at finite 
Pq more reliable. Furthermore, since we have explicit ac¬ 
cess to all fundamental degrees of freedom of QCD, i.e. 
quark, gluon and ghost propagators, we are in a position 
to determine the Polyakov loop potential at all values 
of T and pq and thereby also to study the deconfine¬ 
ment transition. This has been exploited in Refs. [14-16] 
for physical quarks and in Ref. [35] for a study of the 
second-order critical surface of the deconfinement transi¬ 
tion of heavy quarks. These studies are complemented by 
corresponding ones using the functional renormalization 
group, see e.g. [36] and references therein. 


For the gluon DSE, shown in Fig. 1, we use the same 
setup as described in detail in Ref. [14]. We use lat¬ 
tice input for the quenched propagator at different tem¬ 
peratures and determine the temperature and chemical 
potential dependent effects of the quark loop explicitly 
using the quark propagator from its DSE. We work with 
two fermion flavors, Nf = 2, in the isospin limit which 
allows us to use one and the same quark DSE for both 
flavors. The quark loop in the gluon DSE is then sim¬ 
ply multiplied by a factor of two to accommodate for 
both flavors. The essentials of this setup are collected in 
App. B; more details can be found in [14] and shall not 
be repeated here for brevity. 
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FIG. 2. The full, untruncated Schwinger-Dyson equation for the quark-gluon vertex [37] is shown diagrammatically in the first 
equation. The second equation describes the first terms of an expansion in terms of hadronic and non-hadronic contributions 
to the quark-antiquark scattering kernel. In both equations, all internal propagators are fully dressed. Internal dashed lines 
with arrows correspond to ghost propagators, curly lines to gluons and full lines to quark propagators. In the second equation, 
the dotted line describes mesons and the triple line baryons. 
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FIG. 3. DSE for the quark propagator, where the quark-gluon 
vertex is separated into non-baryon and baryon contributions. 


B. Baryon effects in the quark DSE and 
quark-gluon vertex 

Let us now focus on the Dyson-Schwinger equation for 
the quark-gluon vertex and explicate its structure. The 
full equation is shown in the left part of Fig. 2. It con¬ 
tains three one-loop diagrams with a fully dressed quark 
(solid), ghost (dashed) and gluon line (curly) running 
through the loop and attached to the external gluon by 
a corresponding bare vertex. The remaining graph is a 
gluonic two-loop diagram with a bare four-gluon vertex. 
These diagrams contain four- and five-point Green func¬ 
tions that are IPI with respect to the external legs in the 
t channel, i.e., they contain no contributions from inter¬ 
mediate annihilation of the external quarks into a single 
gluon line. The four- and five-point functions can be 
expanded in skeleton diagrams with fully dressed inter¬ 
nal propagators and primitively divergent vertices [37]. 
For our purposes we concentrate on the first non-trivial 
diagram that contains a four-quark amplitude. In its 
skeleton expansion, a part of the resulting diagrams can 
be re-expressed in terms of Bethe-Salpeter vertices and 
propagators of mesons as well as Faddeev-type vertices 
for baryons. 

The result of such an expansion is shown in the sec¬ 
ond equation in Fig. 2. The first diagram corresponds 
to (off-shell) meson exchange between the quark lines. It 
is important to note that this meson is not introduced 
as a new elementary field; it is rather a composite ob¬ 
ject of a quark and an antiquark that is described (at 
least on-shell) by its Bethe-Salpeter equation (BSE). The 
first baryon exchange diagram shows up as a two-loop 
diagram involving the baryons’s Faddeev amplitude. In 
Ref. [38] the corresponding diagram has been displayed 
in the quark-diquark approximation already on the level 
of the vertex DSE; below, we will introduce this approxi¬ 
mation on the level of the quark DSE. Finally, we show a 


representative non-resonant contribution due to dressed 
one-gluon exchange. Note that double-counting is triv¬ 
ially avoided in this combined expansion in elementary 
and effective degrees of freedom due to different quan¬ 
tum numbers in the exchange channel. 

The computation of the hadronic diagrams in Fig. 2 is 
rather involved. The meson-exchange diagram requires 
the solution of a coupled system of the DSE for the 
quark propagator and a corresponding BSE for the me¬ 
son Bethe-Salpeter amplitude. The effect of the pion 
back-reaction onto the quark in the vacuum has already 
been explored to some extent in the context of pion- 
cloud contributions to light mesons and baryons [38- 
40] . Even more complex is the baryon-exchange diagram, 
which involves the computation of a Faddeev-type equa¬ 
tion for the baryon bound state. In the present work 
we are interested in hadronic effects at finite chemical 
potential, which will primarily show up in the diagram 
including baryons because all elements of this diagram 
(the quark and baryon propagators as well as the baryon 
wave function) depend on chemical potential. Since the 
pion-exchange diagram is the hadronic contribution with 
minimal explicit dependence upon chemical potential we 
relegate its explicit study to future work and focus ex¬ 
clusively on the baryonic diagram. 

The resulting quark DSE with an explicit separation 
into non-baryon and baryon parts is shown in Fig. 3. The 
three-loop diagram contained therein is hard to evaluate 
numerically, especially at finite temperature and chemi¬ 
cal potential. To make this diagram tractable, we there¬ 
fore introduce an additional approximation and convert 
the three-quark Faddeev amplitude into a quark-diquark 
Bethe-Salpeter amplitude for the baryon. To this end, 
note that the gluon in the baryon-exchange diagram is 
attached to the quark on the left, but since it couples to 
both quarks symmetrically we can rearrange the diagram 
as illustrated by the first equation of Fig. 4. The incom¬ 
ing and outgoing baryon lines have to be connected by a 
baryon propagator which is indicated by the open circles. 

Inserting a separable quark-diquark ansatz for each 
three-body Faddeev amplitude leads to the second equal¬ 
ity in Fig. 4. The resulting topologies can be grouped 
into two different classes: one where the incoming quark 
couples to a diquark amplitude and one where it cou- 
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FIG. 4. Baryonic diagram in the quark DSE with the quark-diquark approximation of the baryon’s Faddeev amplitude. 
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FIG. 5. Quark DSE with diquark and baryon loop. In these loops the right vertices (circles) are Bethe-Salpeter amplitudes 
whereas the left vertices (hatched circles) are effective ones summing all effects from the diagrams in Fig. 4. In the main text 
we give arguments why these vertices are well approximated by bare ones. 


pies to a quark-diquark amplitude.^ The hatched ampli¬ 
tudes are effective and absorb all the remaining objects 
in these graphs. As a consequence, the quark DSE takes 
the form shown in Fig. 5, which contains a quark-diquark 
and a baryon-diquark loop. For brevity we will refer to 
them as ‘diquark’ and ‘baryon’ loops in what follows. 
In both diagrams the vertex appearing on the right is 
a proper Bethe-Salpeter amplitude, once for a diquark 
and once for a baryon in quark-diquark approximation. 
The hatched vertices on the left carry the same quantum 
numbers as their counterparts on the right but represent 
effective vertices that absorb all effects appearing in the 
multi-loop diagrams in Fig. 4. Thus, the quark DSE in 
Fig. 5 follows directly from the original equation in Fig. 3 
if the quark-diquark ansatz for the Faddeev amplitude is 
made. We will specify its ingredients in Sec. IIC. For 
now, note that we work in the Nf = 2 theory in the 
isospin symmetric limit, which leaves us with isospin sin¬ 
glet scalar and isospin-triplet axial-vector diquarks and 
a degenerate isospin doublet of nucleons. 

For the non-baryonic part of the quark-gluon vertex 
(denoted by ’NB’ in Fig. 3) we employ a construction us- 


^ In principle there are two further diagrams with a closed quark 
loop where the gluon couples to each of the quarks in the loop. 
However, their contributions cancel each other, i.e., the ‘diquark- 
gluon vertex’ is zero. 


ing the first term of the Ball-Chiu vertex that satisfies the 
Abelian Ward-Takahashi identity [41], multiplied with an 
infrared-enhanced function of quark and gluon momenta 
that accounts for the non-Abelian dressing effects and 
the correct ultraviolet running of the vertex. The explicit 
expressions are collected in appendix A. In the following 
subsection we complete the discussion of the quark DSE 
with the last remaining ingredient, the diquark Bethe- 
Salpeter amplitudes together with the quark-diquark am¬ 
plitude for the baryon. 

C. Diquark and baryon amplitudes 

In principle, the baryon is a three-quark state and 
a comprehensive, full treatment of its structure should 
take this explicitly into account. Indeed, the corre¬ 
sponding three-body Faddeev equation has been solved 
in Refs. [40, 42-46] and electromagnetic as well as ax¬ 
ial form factors have been extracted [46-48]. Owing to 
the dynamical formation of diquark correlations inside 
the nucleon, a potentially satisfying approximation to the 
three-body framework is a description in terms of quark 
and diquark degrees of freedom. The BSE for such a 
baryon with quark and diquark constituents is displayed 
in Fig. 6. In this approximation, the quark and diquark 
inside the nucleon interact via quark exchange and the 
corresponding diquark amplitude has to be determined 
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from a separate BSE. 

Using a simple model for the underlying quark and 
diquark propagators and ansdtze for the diquark am¬ 
plitudes, baryon properties in the quark-diquark picture 
have been determined in many works, see e.g. [49-52] 
and references therein. A more fundamental approach is 
the use of an underlying quark-gluon interaction, from 
which all components of such a calculation, the quark 
propagator in the complex momentum plane, the di¬ 
quark amplitude from its BSE, the diquark propagator 
from its scattering equation, and the baryon, are deter¬ 
mined consistently without any introduction of further 
parameters. This has been performed in [53-56] using 
a well-established rainbow-ladder interaction kernel for 
the quark-gluon interaction. In all these calculations it 
turned out that a satisfactory description of the ground- 
state properties of the nucleon and A baryon can be ob¬ 
tained using scalar and axial-vector diquarks only. The 
quark-diquark approximation works well at zero temper¬ 
ature and chemical potential; below we assume that this 
is still the case at finite T and /i^. Whether that is true 
remains to be studied in future work. 

In fact, if one is only interested in the gross proper¬ 
ties of the nucleon even the influence of the axial-vector 
diquark may be omitted and both the diquark and the 
nucleon can be represented by their leading tensor struc¬ 
ture. In this approximation, the diquark and nucleon 
Bethe-Salpeter amplitudes are parametrized by 

P) = fd,{q^) 7 'c 0 ^ , 

^ (3) 

rAr(g, P) = fN{q^) A+(P) 0 ^ 0 C . 

Here, q, P are the relative and total momentum of the 
bound states, C = 7472 is the charge-conjugation ma¬ 
trix, and A-|_ the projection operator onto positive-energy 
states (which we omit in the baryon loop diagram be¬ 
cause its purpose is already served by the nucleon prop¬ 
agator). We use normalized color wave functions with 
capital subscripts and normalized flavor wave functions 
with small subscripts; s° = — du^) = ^* 0-2 with 

Pauli matrices di and = {uu^ + dd^) = 1. The solu¬ 
tions for the diquark and nucleon amplitudes determined 
in the rainbow-ladder framework of Refs. [54-56] are well 
parametrized by 




FIG. 6. The Bethe-Salpeter equation for the baryon in the 
quark-diquark approximation. 

The remaining question concerns the effective ampli¬ 
tudes in Fig. 5 which we did not specify yet. For those 
we resort to a simple approximation: we take them as 
‘bare’, i.e., we use Eq. (3) with fdq{q^) = / 7 v(g^) = 1- 
This is analogous to the treatment of the pion loop in 
the quark DSE in Refs. [38-40]. It can be motivated by 
estimating the overall strength of the baryon diagram in 
Figs. (3-4) from its contribution to the quark conden¬ 
sate. When connecting the quark lines with the scalar qq 
vertex (calculated from its inhomogeneous BSE), the re¬ 
sulting vacuum bubble C is proportional to the integrated 
(off-shell) scalar form factor of the nucleon: 

^ _ 2mN f d^P gs{P'^,Q'^ = 0) 

3 7 (2^)4 P^+m% ■ 

This can be seen by inserting the covariant Faddeev equa¬ 
tion for the three-body amplitude in the first line of 
Fig. 4; the resulting quantity is what appears in baryon 
form factor diagrams such as in Ref. [46]. The on-shell 
value of the scalar form factor at P^ = —m% is de¬ 
termined by the nucleon sigma term via the Feynman- 
Hellmann theorem: 

2 n\ 2 drUN 

dN = mqgs[-mf^,0) = mq- - Ki . ( 6 ) 

dniq dm^ 

The magnitude of C obtained with the experimental 
value dN = 45 MeV, together with an integral cutoff 
aX niN = 0.94 GeV, is similar to the value obtained from 
numerically tracing the sum of the diquark and baryon 
diagrams in Fig. 5 with a scalar vertex, however with 
only one amplitude dressed in each case.^ Dressing both 
would overestimate the strength by far due to the nor¬ 
malization factors in Eq. (4). Since we are only interested 
in the gross effects of baryons on the phase diagram, we 
therefore view this as a justified approximation. 


Pdq \ 


/„(,“) = iV„ («-«■- + 

with X = g^/A^ and the scale A = 0.7 GeV. The normal¬ 
ization factors are obtained from normalizing the corre¬ 
sponding full Bethe-Salpeter amplitudes and are given 
by Ndq = 15.6 and Nn = 28.4. The parameters are 
c^dq = 0.85 and = 1.0 for the exponentials and 
Pdq = 0.02 and Pn = 0.03 for the UV behavior. 


fdqiq^) = Ndq + 


D. Quark DSE including hadronic loops 

Putting everything together, we will now give explicit 
expressions for the diagrams in the quark and gluon DSEs 


^ We also calculated the vacuum bubble from the baryon diagram 
in Fig. 3 directly using a reasonable off-shell ansatz for the three- 
quark Faddeev amplitude instead of a cutoff; the result is in the 
same ballpark. 




















6 


including the hadronic back-reaction diagrams. As will 
become clear below, the hadronic effects on the quark 
propagator enter on the level of ten-percent corrections. 
From the diagrammatic form of the quark-gluon vertex 
DSE this is exactly the order of magnitude as expected, 
since the corresponding diagrams are suppressed by a 
factor 1/A",?. In the Yang-Mills sector of QCD the to¬ 
tal quark effects are on the level of a 1/Nc correction. 
Therefore, hadronic contributions to the quark-loop dia¬ 
gram in the gluon DSE only contribute at 1 /A/? and it is 
well justified to neglect those in a first exploratory cal¬ 
culation. Hence we will use the same truncation for the 
gluon DSE as in previous works [14, 16, 34]. 

In the quark DSE we take into account the three dia¬ 
grams in Fig. 5. If we denote the quark dressing functions 
in Eq. (1) collectively by H(p) = A{p), B (p), C{p) and 
abbreviate the gluon, diquark and baryon-loop contribu¬ 
tions to the quark self-energy by eJ and E^, the 

resulting equations read 

H{p) = Z2 Xh + > (7) 

where Xb = is the bare current-quark mass, = 
Xc = 1, and Z 2 is the quark wave-function renormal¬ 
ization constant. The gluon-dressing loop E|^“® contains 
the unquenched, temperature- and chemical-potential de¬ 
pendent gluon propagator together with a model for the 
quark-gluon vertex [14, 16, 34]; the explicit formulas are 
relegated to App. A. The self-energy contributions from 
the diquark and baryon loop are given by 


^h{p) 


1 fdq{^) Ddqjq + P) dq 

2^ q^A^q)+Q^gCHq)+B^q) 

1 /jv(f -p) Ddq{q-p) 

3 ^ q'^ + {ojq + ‘iipqY + m% ^ 


with 


II 

pz 

T^ba 

II 

"^C{q), 

T^ba 


UJp 


II 

-B{q), 

T^ba 


P Q 

p2 ’ 

UJq + 3ipq 
Cjp 

-tun . 


(9) 


We have already carried out all color and flavor traces. 
At finite temperature and chemical potential the ar¬ 
guments p, q serve as abbreviations for p = (w„,p), 
q = (ujm,q)- The (fermionic) Matsubara frequencies are 
given by LUn = 7rT(2n -I- 1) and we write ui = to + ipq 
with quark chemical potential p,q. The Matsubara sum 
as well as the integration over the loop three-momentum 
q is abbreviated by = I diquark 

propagator Ddq is given by 


Ddq{q±p) 


{q ± pY + {uq ±ujp + 2ipqY + rridg 


■ (10) 


In these expressions we take into account the lowest- 
lying = 1 / 2 + baryon multiplet for the two-flavor 
case, i.e., the nucleon, in the approximation with scalar 
diquarks only. In principle other baryons may also 
contribute but since they are suppressed by powers of 
vaY l'rn?B with respect to the nucleon their influence is cer¬ 
tainly subleading. There is, however, one exception: the 
parity partner of the nucleon becomes (approximately) 
mass-degenerate once chiral symmetry is restored, i.e. 
in the high temperature/density phase. Performing the 
Dirac traces of the corresponding loops, it turns out that 
contributions from mass-degenerate multiplets of parity 
partners cancel each other in E^, while they add up 
in the other two contributions to the quark self-energy. 
We take this effect qualitatively into account by mul¬ 
tiplying the right-hand side of E^ with an additional 
factor M{T, pq)/M{0,0) evaluated at zero momentum 
and lowest Matsubara frequency. Here M = BjA de¬ 
fines the renormalization-point independent quark mass 
function in the medium. This factor has no effect in 
the vacuum but mimics the cancellation of multiplets 
of parity partners in the chirally restored phase, where 
M{T,fiq) becomes small. A corresponding factor of 
2 — M{T, fj,q)/M{0, 0) is added to E/f and E^. Note that 
the diquark loop, which has been derived from the orig¬ 
inal baryon three-body diagram via the quark-diquark 
picture of baryons, contains diquarks only; thus there 
are no damping/enhancement factors in E[^^^ q. 

The remaining unknowns in Eq. (8) are the temper¬ 
ature and chemical-potential dependence of the diquark 
and baryon masses and Bethe-Salpeter amplitudes. Ide¬ 
ally these need to be determined consistently from their 
BSEs evaluated at finite T and pc- This formidable nu¬ 
merical task is yet to be performed and relegated to fu¬ 
ture work. Here, in this exploratory work, we resort to 
the vacuum expressions for the diquark and nucleon am¬ 
plitudes as given in Eq. (4), evaluated at four-momenta 
that include temperature effects in the form of Matsub¬ 
ara frequencies and the results for the masses from the 
corresponding bound state calculations mjsf = 0.938 GeV 
and mdq = 0.810 GeV. Certainly this can only be a first 
approximation on a qualitative level. In order to gauge 
the quantitative effects of including potential changes of 
the baryon and diquark masses and wave functions with 
chemical potential, we will introduce and discuss addi¬ 
tional dependencies on pq in section HI C. 


III. RESULTS 
A. Vacuum 

Before we consider baryon effects at finite tempera¬ 
ture and chemical potential, we first study the impact 
of the different loops on the strength of dynamical chiral 
symmetry breaking in the vacuum. Both diquark and 
baryon loops originate from diagrams in the quark-gluon 
vertex DSE that contain additional quark loops. In gen- 
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FIG. 7. Left: quark mass function with and without diquark and baryon loops included. Right: regularized and normalized 
condensate as a function of temperature with and without diquark and baryon loops. In addition, we display the result with 
rescaled strength in the quark-gluon interaction, see main text for details. 


eral these are chirally restoring, as has been discussed 
in Refs. [38, 57]. On the lattice, the unquenched quark 
propagator has indeed a smaller mass function than the 
quenched one [58]. While this behavior is expected for 
the total sum of all unquenching contributions to the 
quark propagator, it is not necessarily true for each indi¬ 
vidual diagram such as those investigated here. As a first 
exercise we therefore study the sign and magnitude of 
the individual effects of each contribution onto the quark 
mass function. For the calculation we use two dynamical 
quark flavors in the gluon DSE, Nf = 2, a renormalized 
bare quark mass TO(/r^) = 0.8 MeV at the renormaliza¬ 
tion point /i = 80 GeV, and an interaction strength pa¬ 
rameter di = 8.05 GeV^ in the quark-gluon vertex (cf. 
App. A). These values have been taken over from the 
Nf = 2-1-1 theory (set A 2+1 in Ref. [16], matching cor¬ 
responding lattice results) without further adaption for 
reasons of simplicity. 

Our results are displayed in the left diagram of Fig. 7. 
It shows the quark mass function calculated with only the 
usual gluon-dressing loop, compared to results including 
the diquark and baryon loops individually and in combi¬ 
nation. Indeed, all additional contributions are chirally 
restoring, with a larger effect coming from the baryon 
loop. The total contribution of both loops is ^ 15%; 
they reduce the quark mass function at zero momentum 
from M(0) = 640 MeV to M(0) = 545 MeV, whereas the 
impact is immaterial in the large momentum regime.^ 
In total, we find that the effect due to baryons is rather 
large. Lattice QCD finds total unquenching effects in the 


® Note that the considerable size of the generated quark mass is 
a direct effect of performing a Nf = 2-calculation while working 
with scales adapted to the Nf = 2 -|- 1 theory. An additional 
back-coupling of the strange quark would reduce the strength of 
the gluon propagator and decrease the quark mass considerably. 


quark mass function of less then 20 % [58] ; thus our bary- 
onic effects leave almost no room for other unquenching 
corrections like e.g. meson back-coupling effects. This 
may be attributed to the comparably simple approxima¬ 
tion of the baryon wave functions used in this work. On 
the other hand, one of the goals herein is to gauge the sys¬ 
tematic effects of such contributions onto the QCD phase 
diagram. In such a study it seems better to over- than 
to underestimate the induced systematic corrections. 


B. Finite temperature 

Next we assess the effects of the diquark and baryon 
loops on the quark condensate evaluated at finite T. We 
use a regularized expression for the condensate, 

Afh = ( 11 ) 

rrih 

which eliminates the divergences appearing for non-zero 
bare quark masses. The definition of is given 

in Eq. (Cl). Eor the heavy quark mass we choose 
mh{80GeV) = 100 MeV. In order to evaluate the cor¬ 
responding condensate {ip'ip)h in the Nf = 2 theory 
we would need to solve the complete coupled system of 
DSEs, Fig. 1, a second time for each temperature and 
chemical potential. However, the sole purpose of {^ip)h 
is regularization. Thus it turns out to be sufficient to 
evaluate this quantity from the quark DSE with modified 
quark mass mi —>■ in the bare quark propagator 

but keeping the gluon and the quark-gluon vertex (in¬ 
cluding baryonic loops) from the light-quark calculation. 
We have explicitly checked that this procedure is a good 
approximation for some selected values of temperature 
and chemical potential and then adopted it throughout 
the phase diagram. The transition temperatures for the 














chiral crossover are extracted from the maximum of the 
chiral susceptibility. 

Our results for — 0 are displayed in the right dia¬ 
gram of Fig. 7. Compared to the calculation without di¬ 
quark and baryon loops, the additional loops reduce the 
strength of dynamical chiral symmetry breaking and the 
transition temperature for the chiral crossover reduces 
correspondingly. It turns out that this sizable impact of 
the baryonic contributions on the chiral transition can be 
almost completely reabsorbed into the vertex truncation 
by rescaling the strength of the ’NB’-part of the vertex. 
To this end we modify the parameter di (cf. Eq. (A2)) 
such that the critical temperature does not change at 
= 0 upon taking baryon loops explicitly into account. 
This amounts to di = 8.05 GeV^ —>■ di = 8.94 GeV^. 
The resulting condensate (the red solid curve in the plot), 
where baryon and diquark effects are included, recovers 
to very good accuracy the original shape of the conden¬ 
sate. 

This observation is important for our general strategy. 
In Ref. [16] the lattice data for the condensate of the 
Nf = 2 + 1 theory have been reproduced point-wise in a 
formulation using the gluon-dressing loop only, without 
making the baryonic degrees of freedom explicit. Here, 
for Nf = 2, we observe that we can reproduce a similar 
functional dependence of A; ^(T) using explicit baryonic 
degrees of freedom and a rescaled version of the quark- 
gluon interaction. We regard this as a strong indication 
that the same property holds in the Nf = 2 + 1 theory. 
In the following, we therefore use the Nf = 2 theory with 
rescaled interaction strength parameter di as a template 
to study the baryonic effects at hnite /ig. 

While at zero chemical potential all effects can be ab¬ 
sorbed into di, this is not a priori clear for the finite 
chemical potential case because the diquark and baryon 
loops contain a much stronger explicit dependence on pg 
than the gluon dressing loop, as discussed above. In the 
next section we will explore the consequences of these ad¬ 
ditional contributions for the location of the critical end 
point (CEP). 


C. Finite temperature and chemical potential 

In the Dyson-Schwinger approach to the QCD phase 
diagram the introduction of (real) quark chemical poten¬ 
tial is straightforward; cf. Eq. (I), where the dressing 
functions of the quark propagator become complex. The 
impact of chemical potential is apparent in the quark 
propagator but it also affects the gluon (cf. App. B) as 
well as the quark-gluon vertex (cf. App. A) due to the 
explicit unquenching procedure. In this section we in¬ 
vestigate how this nontrivial influence is modified by the 
baryon and diquark loops. We use the maximum of the 
chiral susceptibility, Eq. (C3), as the definition of the 
(pseudo-) critical temperature. 

In Fig. 8 we show the results in the T-pg plane. The 
solid (black) curve is the result for the unquenched sys- 



FIG. 8. Comparison of the phase diagram for Nf=2 including 
different types of selfenergy contributions. 

tern with gluon-dressing loop and no baryonic effects. For 
small values of the chemical potential pg = 0 the transi¬ 
tion is a cross-over up to the filled circle, which indicates 
the 2nd order critical end point (CEP) at the critical 
value pq = Pg. The two dashed lines emerging from the 
CEP mark the first order spinodal region for pg > Pg. In 
this case we find a critical endpoint at 

(r=, ^=) = (177,134) MeV (12) 

for the two-flavor theory Nf = 2. Comparison with the 
dashed (red) curve, which includes explicit baryonic ef¬ 
fects with rescaled vertex strength, leads us to the fol¬ 
lowing observations: 

• a critical endpoint still exists; 

• the chiral phase transition lines are almost on top 
of each other; 

• the critical endpoint is shifted by less than 5 MeV 
to smaller chemical potential. 

The first observation is important because in QC 2 D a 
disappearance of the critical endpoint was observed af¬ 
ter introducing the two-color equivalent of baryonic ef¬ 
fects [27, 28], see also [29-31]. While a negligible in¬ 
fluence of baryonic degrees of freedom on the transition 
line is generally expected at low chemical potential, in 
our case such a behavior also persists for higher pg and 
remains true for the critical endpoint. This also implies 
that our original truncation with the unquenched gluon¬ 
dressing loop can implicitly absorb baryonic effects at 
finite chemical potential, at least those that are captured 
by our simple approximation. 

Potentially important effects beyond our current 
scheme are additional dependencies of the baryon’s mass 
and wave function on the chemical potential. Unfortu¬ 
nately not much is know in this respect. In Ref. [60] the 
authors investigate the thermal properties of baryons at 
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/ig = 0 and find that the amplitudes are almost indepen¬ 
dent of T, whereas the masses rapidly increase around 
the (pseudo-) critical temperature. As a result, baryonic 
contributions would decrease in importance in a region 
close to the transition line. Such temperature-dependent 
effects would presumably have no impact on our results 
since they can be reabsorbed in the strength of the ’NB’- 
part of the vertex. 

In the absence of explicit knowledge, we gauge the im¬ 
pact of modifications of the baryon wave function with 
chemical potential by multiplying the baryon loop with 
a function 


= 1 - 




(13) 


where we make use of recent evaluations of the curvature 
of the chiral transition line on the lattice [9-11]. This 
curvature can be parametrized in terms of a quantity k, 


T(0) 


1 — AC 


T{pb) 


2 


(14) 


that characterizes the lowest order in a Taylor expan¬ 
sion in the baryonic chemical potential While recent 
lattice values for k range between 0.0135... 0.020, our 
value is somewhat larger. Neglecting 1/A^/ corrections, 
we adopt the Nf = 2-1-1 value n = 0.0149 from [10] 
for our Nf = 2-calculation and match the coefficients 
AK,a„,6„ in Eq. (13) such that we reproduce the lat¬ 
tice curvature in a region where the lattice can well be 
trusted. This is possible for = 0.714 GeV, = —10.3 
and 6 k = 36. The function /«, then has a minimum 
at fig « 120 MeV. For larger chemical potential we use 
fnip-q) = = 120MeV) to make the function mono¬ 

tonic. The resulting phase diagram is shown as the dash- 
dotted (indigo) curve in Fig. 8, with a new location of the 
critical end point at 


(T^/r=) = (197,125)MeV. (15) 

Due to the smaller curvature at low chemical potential, 
the CEP shifts by ~ 10% towards larger temperature 
and ~ 5% towards smaller chemical potential. The ra¬ 
tio changes accordingly from = 2.3 to 

= 1.9. These changes are by no means dramatic 
but quantitatively significant. However, this comes at the 
expense of a modification of the strength of the baryon 
loop by more than 50% due to the additional function 
/„. Whether such a variation of the baryon wave func¬ 
tion and masses with chemical potential is realistic or not 
needs to be investigated in the future. It is also by no 
means clear whether baryonic effects are the only possible 
source of a smaller curvature in the present framework. 
Here we only demonstrated that it is possible in principle 
that baryonic corrections can induce such an effect. 

Qualitatively, it is always the case that the CEP shifts 
towards larger temperatures and smaller chemical poten¬ 
tial if fnifJ-q) is smaller than one for > 0. The opposite 


effect can be obtained if fKifJ^q) is chosen to be larger than 
one: the CEP then shifts towards smaller temperatures 
and larger chemical potential. Apart from arguments by 
comparison with the lattice we see no physical reason 
a priori why baryon effects should have one effect or the 
other. Again, this needs to be studied in a more advanced 
framework. 


IV. SUMMARY AND CONCLUSIONS 


In this exploratory study we extended our existing 
truncation of the coupled system for the quark and gluon 
Dyson-Schwinger equations to take explicit baryonic de¬ 
grees of freedom into account. This was achieved by 
considering a specific class of diagrams in the Dyson- 
Schwinger equation for the quark-gluon vertex where 
genuine hadronic contributions can be identified. Upon 
introducing baryons through the quark-diquark picture, 
the baryon diagrams enter as a baryon-diquark loop and a 
quark-diquark loop in the quark Dyson-Schwinger equa¬ 
tion in addition to the gluon-dressing loop. In the Nf = 2 
calculation performed herein we employ the vacuum am¬ 
plitudes and masses for the nucleon and the (scalar) di¬ 
quark, but we take into account cancellation effects due 
to the degeneration of the chiral partner of the nucleon 
through a factor that couples the baryon-diquark loop 
to the chiral dynamics of the system. With this setup 
we performed a calculation of the QCD phase diagram 
and find that the inclusion of baryon degrees of freedom 
changes the location of the critical endpoint only by a 
few MeV in T and ^,q. 

More drastic effects are possible once a dependence of 
the baryon masses and wave functions on chemical po¬ 
tential are taken into account. We estimated these using 
a parametrization that reproduces the lattice transition 
line at small chemical potential. As a result we find a shift 
of the CEP in the 5 ... 10% range which drives the ratio 
slightly below 2. We expect that these results ob¬ 
tained in the two-flavor theory still hold qualitatively for 
the Nf = 2-1-1 case. This will be explored in future work, 
where we also strive to determine the chemical potential 
dependence of the baryon and diquark masses and wave 
functions explicitly. 
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Appendix A: Gluon contribution to quark 
self-energy 


and the kernels K read 


To make the paper self-contained we collect in this ap¬ 
pendix the ingredients of the non-baryonic (NB) part of 
the quark DSE, i.e., the gluon-dressing loop in Fig. 5. 
Our ansatz for the quark-gluon vertex that appears 
therein is given by 


g; k) = 7 ^ Ts -t 5^4 r4) r(/c2), (Al) 


where 


Ts 


r(fc 2 ) 


A{p) -b A{q) 

" 2 

di X 

d2 + k^ cc -b 1 


r4 = 
l3oa 

47r 


C{p) + Ciq) 


ln(l -b x) 


2S 


(A2) 


Kaa 


Kac 

Kbb 

Kca 

Kcc 


hs 


f ■ k q ■ k 




+ 2 


(rs 


fc 2 A :2 

Zt p ■ k q ■ k 


fc 2 fc 2 

, Zt 

'■p" 


Zj' 

“p" 


fc 2 


r 

fc 2 /fc 2 P ■ 


A:2 A:2 


Zl (jjI 

^ -b r4 


Zl q k 

Wfe , 

k^ fc2 

, Zl ujI 


^ k^ fc2 

jUJq- 


P * 2 ^ 


/c 2 /fc 2 


(A4) 


Appendix B: Unquenching the gluon 


with a; = fc^/A^. Here, p = {ujp,p) and q = {u}q,q) 
are the fermion momenta and k = {(jjk,k) is the gluon 
momentum. The scales ^2 = 0.5 GeV^ and A = 1.4 GeV 
and the coupling a = 0.3 are adapted to match the corre¬ 
sponding scales of the quenched lattice gluon propagator 
that we use in the gluon DSE. Whereas ^2 and A control 
the renormalization-group running of the vertex function 
from the large- into the low-momentum region, di con¬ 
trols the strength of the quark-gluon interaction at small 
momenta and therefore the amount of dynamical chiral 
symmetry breaking in the hadronic phase. Its value is 
discussed in the main text. The ultraviolet momentum 
region is governed by the anomalous dimension of the 
vertex <5 = -9-, 


iiN^-SNf ^“^d /3o — 


llAr„-2Afj- 


The explicit expressions for the gluonic parts of the self 
energy are: 


y^glue 


■^glue 


y^glue _ 

- 


Z2CFg^^ 

Q 

Z2CFg^^ 

9 

Z2CFg^^ 


r(fc") A{q)KAA + C{q)KAC 

D{q) p2 

-jj^B{q)KBB, (A3) 

r(fc^) A{q)KcA + C{q)Kcc 

D{q) Qp 


In the Yang-Mills sector we solve the gluon DSE in¬ 
cluding the quark-loop contributions shown in Fig. 1, i.e. 


Nf 

D-Ak) = [D^Zik)]-^ . (Bl) 

/ 

^iAk) = ^ Tr [ 7 ^ q-, k) 5 ^( 5 )] , 


with the explicit flavor dependence indicated by the su¬ 
perscript /. DZ^-{k) denotes the quenched gluon propa¬ 
gator which is taken from lattice calculations [59] ; the fit¬ 
ting procedure has been discussed in Ref. [33]. r^(p, q\ k) 
is given in Eq. (Al) but T is evaluated for p^ -b q^ 
instead of Z to ensure multiplicative renormalizability 
(cf. Eq. (B4) below). The lattice fits for the gluon dress¬ 
ing functions Zt.l (cf. Eq. (1)) are given by 


ZT,L{k) = 


[x + iy 


~ \bT,LiT) 

X + aT,LiT)) 

^ln(l+x)'" 


(B2) 


where k = p—q, Cf = | is the Casimir operator, Z 2 is the 
wave-function renormalization constant, and T has been 
defined above. The Matsubara sum as well as the inte¬ 
gration over the loop three-momentum q is represented 
by = T Xn, I (Arfz' "bhc denominator of the quark 
propagator is given by D(g) = q'^A'^{q)+cOq C'^{q) + (q) 


where x = A /and aT,L(T), hT,L{T) are temperature- 
dependent fit parameters. In the ultraviolet, the log¬ 
arithmic term leads to the perturbative running with 
anomalous dimension 7 = 22 n'"-anI ■ "bbe temperature- 
independent parameters are c = 5.87 and A = 1.4 GeV. 
With a transition temperature of = 277 MeV for the 
quenched SU(3) theory, the temperature-dependent pa- 
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rameters are given by 


aL{t) 


arit) 

bhit) 

brit) 


0.595 - 0.9025 • t + 0.4005 if i < 1 

3.6199 • t - 3.4835 if f > 1 ’ 

0.595 + 1.1010-^2 if t < 1 

0.8505 • t - 0.2965 if f > 1 ’ 


1.355 - 0.5741 • t + 0.3287 iit<l 
0.1131-^ + 0.9319 if f > 1 ’ 


1.355 + 0.5548-^2 if t < 1 

0.4296 -1 + 0.7103 if f > 1 


(B3) 


with t := T/Tc- Note that since this expression repre¬ 
sents the quenched gluon propagator it is independent 
of chemical potential, which enters the gluon DSE only 
through the quark loop, and Nf = 0 in the anomalous 
dimension. 

By contracting Eq. (Bl) with the projectors in Eq. (2) 
one arrives at the equations for the transverse and longi¬ 
tudinal parts of the quark loop. Since we are using hard 
momentum cutoffs in the numerical integration, these 
need to be carefully regularized to remove quadratic di¬ 
vergencies without spoiling the Debye screening masses. 
This procedure is described in the appendix of Ref. [14] . 
The resulting equations read 


(B4) 


nT,L(fc^ +^fc;0)- We use the same approximation in 
the quenched gluon propagator where we only have the 
zero mode from the lattice. This procedure is surpris¬ 
ingly accurate, as has been verified e.g. in [33, 61]. 


Appendix C: Chiral condensate and definition of Tc 
in a crossover region 


The quark chiral condensate for a flavor / is given by 

W)/ = j ^Tr,5 [Sfip)] , (Cl) 


where Z 2 is the quark wave function renormalization 
constant, Z^. the quark mass renormalization constant, 
Nc the number of colors and Trp indicates the Dirac 
trace. Due to the nonzero quark mass the condensate is 
quadratically divergent and has to be regularized. We 
consider two ways to define a (pseudo-) critical temper¬ 
ature that are both connected to the chiral condensate. 
The first one is the inflection point of the chiral conden¬ 
sate with respect to temperature. 


rpinfl. 


= max 
VT 


dT 


(C2) 


for the contribution with lowest Matsubara frequency, 

where p and q = p + k are the quark momenta in the and the second one returns the maximum of the chiral 
loop and susceptibility: 


KT = A{p)A{q)rs{3 


Kl = A{jp) A{q) 


{p ■ kf 


fc 2 


2p ■ k — p‘ 


, pkkq 

(2-p- p q 


Tip - q 


+ B{p) B{q) (r4 - Ts) - C{p) C{q) (Ts + r4) c^l. 


The higher modes are accessed via Bxxik^ 


= max 


dnrif 


(C3) 


Both definitions are independent of the regularization of 
the divergent part (finite bare masses). While they give 
different results for Tq in a crossover region, they return 
the same position for the critical end point. 
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